Non-perturbative approach to the quantum Hall bilayer 
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We develop a non-perturbative approach to the quantum Hall bilayer (QHB) at = 1 using trial 
wave functions. We predict phases of the QHB for arbitrary distance d and, our approach, in a 
dual picture, naturally introduces a new kind of quasiparticles - neutral fermions. Neutral fermion 
is a composite of two merons of the same vorticity and opposite charge. For small d (i.e. in the 
superfluid phase), neutral fermions appear as dipoles. At larger d dipoles dissociate into the phase 
of the two decoupled Fermi-liquid-like states. This scenario is relevant for the experimental situation 
where impurities lock charged merons. In a translation invariant (clean) system, continuous creation 
and annihilation of meron-antimeron pairs evolves the QHB toward a paired phase. The quantum 
fluctuations fix the form of the pairing function to g{z) = 1/z* . A part of the description of the 
paired phase is the 2D superconductor i.e. BF Chern-Simons theory. The paired phase is not very 
distinct from the superfluid phase. 

PACS numbers: 73.43.-f, 73.43.Nq, 03.65.Vf 
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I. INTRODUCTION 

The quantum Hall bilayer at i> ~ I consists of two 
layers of two-dimensional electron gases that are brought 
close to one another in the quantum Hall regime of strong 
magnetic fields. When the distance between the layers is 
much smaller than the average distance between electrons 
inside each layer, inter and intra Coulomb interactions 
are about the same. Then the expected v = 1 state is 
the state of a single layer filled lowest Landau level (LLL) 
generalized to two species. There is obvious degeneracy 
in dividing electrons into two groups which leads to the 
phenomenon of spontaneous symmetry breaking [l| and 
the existence of a Goldstone mode [!]■ The expected su- 
perfluid behavior was verified also by very large zero bias 
voltage peak in tunneling conductance [3|, but no clear 
evidence was found for finite temperature Berezinskii- 
Kosterlitz-Thouless (BKT) transition [3| in transport ex- 
periments [5;]. 

Therefore there is a need to systematically address the 
question of superfluid disordering in the QHB. In partic- 
ular there is a need to understand the role of quantum 
disordering in this system that becomes important as the 
distance between the layers is increased. In most of the 
previous work the starting point for the discussion of the 
physics of the bilayer was the ground state (GS) for very 
small distance between the layers as a mean field solution 
to which none or some corrections were developed [1, @ . 
We will take a non-perturbative approach inspired by the 
Laughlin solution of the i/ = 1/3 problem in which we 
will uniquely determine possible wave functions (WFs) 
for the GSs of the bilayer at an arbitrary distance. 

There are two basic paradigms of superfluid disor- 
dering that are known: (1) BKT (2D XY model) for 
which the transition proceeds via unbinding of dipoles 
of vortex-antivortex pairs, and (2) A transition type (3D 
XY model) for which the transition is characterized by a 
condensation of vortex-antivortex loops 0] ■ 



On the other hand, in this paper, through an analysis 
of the allowed possibilities for homogeneous WFs as the 
distance is varied, we will identify two families of WFs 
and relate them to the two ways of disordering the QHB 
superfluid mentioned previously. The families will be in- 
troduced in Section II. 

One family, as it will turn out does not include el- 
ementary vortices - merons of QHB, in its description 
of superfluid disordering. Merons are part of the de- 
scription of the QHB superfluid for small distances as 
is well-known and well-established in Ref. [1]. There- 
fore this family of (homogeneous) WFs can be relevant 
only for dirty systems - systems with impurities, which 
can lock merons due to merons being charged quasipar- 
ticles. Then the only vortices that may participate in 
superfluid disordering and on which description of this 
family of WFs is based, are neutral composites of two 
merons of opposite charges - neutral vortices, and as we 
will find fermionic quasiparticles that carry only layer 
degree of freedom. We will show that the superfluid 
disordering of this family can be understood through a 
Coulomb (fermionic) plasma picture of dipoles of these 
neutral fermions. Therefore this family we can consider 
as the one that exemplifies the BKT way of superfiuid 
disordering, our first paradigm. This whole picture will 
be corroborated by the fact that the WFs of this fam- 
ily do not incorporate quantum fiuctuations (Section IV) 
and, therefore, do not incorporate quantum disordering 
that is based on merons. The family from the viewpoint 
of a dual description (i.e. in terms of quasiparticles - 
neutral fermions) will be analyzed in Section HI. 

The other family incorporates weak pairing among 
neutral fermions and, as we will show, by assuming a 
special kind of pairing agrees and correlates with the de- 
scription of quantum fluctuations of the usual superfluid 
disordering in a translatory invariant system that one 
finds in other approaches (field-theoretical). It is ex- 
pected that this kind of disordering and pairing would 
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lead to a charge density wave (CDW) solution Still 
our general considerations open possibilities for other 
kinds of weak pairing that can be present in this quan- 
tum Hall system. The most likely candidate is the one 
with pairing function g{z) ^ ^ that results in non-trivial 
corrections (from quantum fluctuations and disordering) 
to the ground state wave function as the distance is var- 
ied. In general we expect that a weak pairing scenario 
will correspond to the superfluid disordering of the usual 
superfluid in 2 -I- 1 dimension and therefore to the class 
of 3D XY, our second paradigm. The family with weak 
pairing ansatz will be analyzed in Section IV. 

With respect to the experiments, where impurities are 
necessarily present and we can expect also inhomoge- 
neous ground state solutions, our homogeneous candi- 
dates of the first family (without neutral fcrmion weak 
pairing) are still possible solutions for which transitions 
may proceed via dissociation of dipoles - pairs of oppo- 
site vorticity neutral fermions. In this sense and as will 
be more clear later, the quantum phase transitions with 
respect to changing the distance in Refs. [l,[l3,[lH corre- 
spond to this dissociation. On the other hand, an analy- 
sis will show that in a translatory invariant system meron 
excitations via their loop condensation may produce an 
intercorrelated, paired liquid state for the neutral sector, 
if a transition to a CDW does not occur. 



II. UNIVERSALITY CLASSES OF GROUND 
STATES 

A. Introduction 




(a) (b) 




FIG. 1: Universality classes of wave functions 



filling factor 1/2 WFs; each describes a Fermi-liquid-like 
state flE^ . 

*i/2(H = V{Ts{w,w) - "'^t)'} (2) 

where J^s is the Slater determinant of free waves of nonin- 
teracting particles in zero magnetic field and V represents 
projection to LLL. Underlying quasiparticles are compos- 
ite fermions (CFs) , the usual quasiparticles of the single 
layer quantum Hall physics. 



B. Two families - universality classes of wave 
functions 



A great deal is known from the experimental and theo- 
retical point of view of the QHB in the two extremes when 
the distance between layers, d, is (1) much smaller or (2) 
much larger than the magnetic length, Ib = {h/eBy^^, 
where B is - the magnetic field, the characteristic dis- 
tance between the electrons inside any of the layers. 
When d << Ib, i-e. inter and intra Coulomb interac- 
tions are about the same, the good starting point and 
description is so-called (111) state fl^ . 

i<j k<l p,q 

(1) 

where Zi^ and are two-dimensional complex coordi- 
nates of electrons in upper and lower layer respectively 
and we omitted the Gaussian factors. This is suggestive 
of the exciton binding [isj ; any electron coordinate is also 
zero of the WF for any other electron coordinate - the 
correlation hole is just opposite to electron. This exci- 
ton description can be a viewpoint of the phenomenon of 
superfiuidity found in these systems @, ll| and is closely 
connected to the concept of composite bosons (CBs) [lj| 
that can be used as natural quantum Hall quasiparticles 
in this system. When d >> Ib have the case of the 
decoupled layers and the GS is a product of single layer 



To answer the question of intermediate distances we 
may try to, classically speaking, divide electrons into two 
groups, one in which electrons correlate as CBs and the 
other as CFs [l^. The ratio between the numbers of 
CBs and CFs would be determined by the distance be- 
tween layers. The WF constructed in this way would 
need an overall antisymmetrization in the end, but also 
intercorrelations among the groups as each electron of the 
system sees the same number of flux quanta through the 
system (equal to the number of electrons). This requires 
that the highest power of any electron coordinate is the 
same as the number of electrons in the thermodynamic 
limit. If we denote by a line the Laughlin-Jastrow factor 
Y[j^ b(-^a — zb) between two groups of electrons, A and B 
{A, B = CB, CF), the possibilities for the QHB GSWFs 
can be summarized as in Fig. 1. 

If we ignore the possibility of pairing between CFs [l7j 
denoted by wriggly lines in Fig.l(c)and Fig. 1(d) we have 
two basic families of the GSWFs depicted in Fig. 1(a) and 
Fig.l (b). The requirement that each electron sees the 
same number of flux quanta through the system equal 
to the number of electrons (we are at v — 1) very much 
reduces the number of possible states - wave functions 
in the mixed CB - CF approach. We can consider, for 
example, the possibility (a) depicted in Fig. 1 which 
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stands for the following wave- function in the LLL, 

^-1 = VA1Al{]J{z,^ - zjt) Y[{zki - zii) Y[izp^ - zgi) 

i<j k<l p,q 

i<j 
k<l 

i,j k,l 

where Ai and Ai denote the overall antisymmetrizations. 
In the thermodynamic limit, the relation between the 
number of particles and flux quanta reads, 

iV^ = m^+Nbi + Nf^ + Nf<^, 

= 27V/i+7V,t+^H, (4) 

where we denoted by and N^" separately the number 
of flux quanta that electrons that correlate as CBs and 
CFs respectively see, N^a and Nf„ are the number of 
CBs and CFs inside the layer cr, respectively (cr =ti i is 
the layer index) . The requirement constrains — = 

N^"' , where is the number of flux quanta through the 
system. (This leads to the additional requirement iV/| = 
Nfi which leaves iVb^ — Nti unconstrained, connected 
with the Bose condensation phenomenon that the wave 
function should be part of ( [j, [H, ) . 

The only additional way to count the flux quanta that 
electrons see, with the (symmetric under t"^! reversal) 
application of the Jastrow - Laughlin factors that we need 
to have, is: 

iVf = 2iV/x + 27Vfci, (5) 

which leads to the possibility (6) (with constraints iVf,| = 
Nbi and Nf\ — Nfi) The intercorrelations in the flrst 
family in Fig. 1(a) are in the spirit of ^'m correlations, 
and those in the second family in Fig. 1(b) are in the spirit 
of the decoupled state, '^1/2 x ^'1/2, where we correlate 
exclusively inside each layer. 

C. Discussion 

We can imagine a mixture of both intercorrelations (of 
Fig. 1(a) and Fig. 1(b)) in a single wave function but 
these mixed states, by their basic response [l^, fall into 



one of the universality classes depicted in Fig. 1. In Ref. 
[Tsl explicitly such a mixture and possibility under name 
"generalized vortex metal" was considered, in the scope 
of a Chern - Simons (CS) theory, and it was proved that 
it does not support a Goldstone (gapless) mode which 
was found to exist for the state depicted in Fig. 1(a). 
These generalized states belong to the universality class 
of the state depicted in Fig. 1(b) for which in the scope 
of the same theory we find in the low-energy spectrum 
only a gapped collective mode [ll| . 

The Chern-Simons theory we mentioned neglects the 
overall antisymmetrization built in the classes of Fig. 1. 
We can justify this neglect (1) by taking a point of view 
that stems from similar situations with quantum Hall 
states like hierarchy and Jain's constructions that in the 
low-ener gy sector can be considered as multicomponent 
systems [1^ (we will argue later that the state of Eq. ([3]) 
can be mapped to a hierarchy construction), or (2) a 
posteriori because the results of the effective description 
of the classes in Fig. 1 are quite sensible and are ex- 
pected for the states we are familiar with from numerics 
(the state in our Fig. 1(a) as analyzed in [11]). (We 
do not ask this type of theory for detailed answers any- 
way.) In this way it was found by us (RefslH and [20) 
examining the basic response in the pseudospin chan- 
nel in the random phase approximation (RPA) of these 
Chern - Simons theories that the states in Fig. 1(a) and 
Fig. 1(c) represent superfluids, and the states in Fig. 1(b) 
and Fig. 1(d) represent disordered superfluids, compress- 
ible and incompressible, respectively. (Later, in a more 
complete study, we will find that the states of Fig. 1(d) 
are also compressible in the neutral channel.) 

The two basic possibilities of connecting two extremes 
as depicted in Fig. 1, i.e. without and with pairing 
of CFs, must correspond to the two possible ways or 
paradigms that we know of disordering a superfluid. We 
will substantiate this claim further by examining the two 
superfluid constructions (Fig. 1(a) and Fig. 1(c)) in more 
detail. 



III. NEUTRAL FERMIONS AND BKT 
DISORDERING 

A. Dual picture of the first family of wave 
functions with neutral fermions 

Let us write out the unprojected in the LLL version 
of the construction in Fig. l(a)(Eq.(I3])) in the following 
way: 

■^i=A^Ai { 1'iii(zT.^i)*i/2(u'T)*i/2K) 
X n^^^T - Wjt) Wi^k] ~ wii) 

i,j k.l 

X Yliz^l - Wgt) YliZnil - Wnl)}, (6) 
p,q m,n 
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where, as before, ^^'s and Wa-'s denote coordinates of 
electrons belonging to the layer with index a and and 
Ai, as before, stand for the antisymmetrizations. Us- 
ing 5| and 5j, symmetrizers inside each layer, the same 
function, ^'i, can be written as: 

e C (llk<liWk^ -'Wll)I\p<qiWpl-Wql) 



Our first question may be why states as 



(10) 



would not make a bosonic basis. We look for the following 
overlap: 



111 



(7) 



In the expansion of the density sums we may get 



(11) 



where ^'m denotes the Vandermonde determinant 
(Slater determinant in the LLL) of all coordinates in 
which all groups equally participate. 

By using the expressions for the densities of electrons 
in each layer, p'^(ry) = J2i ^'^iv ~ -^f )j here now Zo- 's de- 
note all electrons of the layer a, we can rewrite the wave 
function in the following way. 



^-1 



d rii] 



Uk<liVk-{ - mi) Up<qiVpl - Vql) 



6'{V[ - zl)6'{v'2 ~ 4)S\vi - 4)S\m - 4) • • • (12) 

which would lead to the following contribution after z- 
integration: 

s'irA - v'2)^'ivi - - m\' + !%!')} 

(13) 



\vi - m\ 



■exp{-(|?7i| 



pT(r/i|)---pJ-(r;„J*iii(z|,Z|) 



where n is the total number of electrons that correlate 
as CFs. The equality is exact; any time we have in the 
product of p's the same layer electron coordinate more 
than once, the Laughlin - Jastrow factors of ?7's in the 
same layer force the wave function to become zero. The 
expression in Eq. ([5]) reminds us of a dual picture in terms 
of some quasiparticles with rj coordinates as in j2l| . Cer- 
tainly we are not describing the incompressible physics 
of a Laughlin state where quasihole operators of coherent 
states span the basis of low-energy physics and allow the 
description in terms of wave functions of quasiholes (dual 
picture) [2lj . Nevertheless we will argue that we can de- 
lineate a sector (find a subspace) to which constructions 
(Eq.® where n is arbitrary) belong, which is spanned 
by a quasiparticle basis of some neutral fcrmionic quasi- 
particles. 

To find those quasiparticles we will rewrite Eq.® as 



The last term, before the dots, comes after the integra- 
tion over z's that do not participate in the delta func- 
tions. As usual [2l[ the term is the result of the screen- 
ing of plasma which we find in the plasma analogy of 
(S^fiii state in its charge channel. The term exactly can- 
cels the preceding one (it is equal to its inverse) and the 
same cancelation will happen for any pair of 77's (in the 
place of • • •) that in remaining z- integration have the role 
of impurities (of charge one) in the plasma of remaining 
z's. This is very good because of our goal to find basis 
states and leaves us to consider only delta functions in 
the contribution. But we can see immediately in Eq. (|13|) 
that 5{r]'-^ — 773) spoils our goal that the states mimic a 
Fock basis of bosonic particles. Therefore as candidates 
for basis states we should consider fermionic states: 

1 



hlT ■■■^nl >= I exp{i(/)(?7i| • • • X 




^-1 



d'^rm 



d^n. 



|exp{i(/)(77i|---?7„j^)} X 
pHvh)- ■ ■ pKvni)'^iii{zi,zi 



where exp{i(j){r])} factor denotes the phase part of the 
Laughlin-Jastrow factors in front of the Fermi seas in 
Eq. ([5]) . With respect to Eq. ([5]) we are allowed to take for 
definiteness that the phase factor always vanishes when 
any of two jy's (or more) from the same layer coincide. 



pH'7it)---pH^M)I*iii >(14) 

for which we can not get contributions of the type in 
Eg. p^ because the phase part does not allow two (or 
more) quasiparticles to coincide. ( Eq. p4)) represents a 
IMJfermionic state for rj quasiparticles because of the phase 
part introduced in Eq.(9) which is antisymmetric un- 
der the exchange of 77's.) Therefore we should con- 
sider fermionic states in Eq. (|14p because of the previ- 
ously found non-desirable terms in the bosonic case (we 
are looking for quasiparticles and their basis states that 
would have features of the Fock space basis): the terms 
(gjike the one with S{ri-^ — %) in Eq. (fT5)) lead to the ab- 
sence of orthogonality of these states which we would like 
to represent coordinate basis states in the bosonic case 
and that can be mended by taking fermions - then these 
terms are absent. By a similar analysis which lead to 
Eq. p3p . considering various possibilities for delta func- 
tion contributions of density operators we can find that 
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the leading - most singular and coherent behavior of the 
states defined in Eq. fH]) is 



< '7lT>%T ■ ■■Vnl\mhV2] ■ --Vnl >- 



(15) 



The rest of contribution constitute incoherent phase fac- 
tors with fewer number (< n) of delta functions but of 
the same kind as in the leading behavior. We can not 
prove that the states make exactly a Fock space of neu- 
tral fermionic quasiparticles, i.e. we do not have an ex- 
act equality in Eg. p^ . but they stand fairly close to 
that status. In other words we do not have the exact 
equality in Eq.(15) i.e. equality to delta functions only 
(appropriately antisymmetrized) , but we have in addi- 
tion some finite contributions which can not change the 
fact that the overlap is singular - at its maximum when 
T] 's and r?'s coincide. Therefore quasiparticles are not 
point-like fermionic quasiparticles (one certainly can not 
expect that from quasiparticles in a strongly correlated 
system), they are extended, but clearly the overlap has 
the singular contribution of antisymmetrized delta func- 
tions which points out that we are fairly (to a good ex- 
tent) close to the fermionic Fock basis description. Even 
in the Laughlin case we can not prove the exact LLL 
delta function overlaps of coherent states of quasiholes. 
The quasiparticles are neutral because in the construc- 
tion of the states there is no net magnetic flux through 
the system. See Ea.([T^ and the definition of the phase 
factor in Eq.® with Eq.®. 

Now that we know basis states, just by looking at 
Eq.® we can read out the GSWF in the dual picture 
in terms of neutral fermions. 



^diial{ri) 



Y\k<i \Vk\ - m\\]\p<q \ripi - riqi\ 



(16) 

This is a wave function of a 2D Coulomb fermionic 
plasma. In the literature 2D Coulomb fermionic plasma 
with same charge particles is fairly known and explored 
[12, • It is a dynamical system of fermionic particles in 
2D that interact with the long-range (^ — ln{r}) interac- 
tion. As shown in Ref. [13 the Jastrow factor of the type 
T\i<j (7 proportional to the interaction coupling 

constant), together with multiplying Slater determinant 
of free waves, describes the ground state function in the 
long-distance limit. In our case we have a generalization 
of such a system to the one with opposite charges. As- 
suming that the concentration of particles is not large, 
which is the case of interest to us, we then expect the 
dipole configurations of particles that the wave function 
in Eq. (fTB|) describes. 
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FIG. 2: The quadrupolar configurations of merons that make 
neutral fermion pair. Compare the same configuration of 
Laughlin quasiparticles as a description of "magnetophonon" 
branch in the Laughlin case in Ref. [T^ . 



B. Discussion 



Merons are true elementary vorticity quasiparticles of 
the translatory invariant QHB system at least for small 
distances between layers as shown in Ref. 4] and carry 
both charge and vorticity. Therefore the neutral fermion 
basis that we described can be a complete basis for the 
ground state evolution of the QHB in the non-translatory 
invariant case in which merons by their charges are bound 
to impurities. 

The wave function in Eg. pop describes the superfluid 
state in Fig 1(a). It encodes dipole positioning of op- 
posite vorticity (layer index) neutral fermions. With 
increasing distance there are more dipoles of neutral 
fermions and they are expected to be less tightly bound 
as in the description of a BKT disordering of a 2D system 
with increasing temperature. Therefore we do not find 
quantum fluctuations in this case. This will be explicitly 
shown by calculations in the following section (see also 
Appendix A). 

In the superfluid phase, with respect to merons, a neu- 
tral fermion dipole should be in essence a superposition of 
quadrupolar combinations of merons - two dipoles which 
come in pairs but at arbitrary distance as illustrated in 
Fig. 2. In this way, as special configurations of dipoles, 
neutral fermions, we expect, constitute the lowest lying 
states of the QHB - (pseudo)spin or phonon waves [^.Tl3|. 
If neutral fermions may be considered as eigenstates they 
must lie very high in spectrum; like electrons in fractional 
quantum Hall states they constitute the physics of ^'i but 
their wave function Eq. (fT6|) describes a highly correlated 
state. 

The dual expression of Eq. p6p was derived under as- 
sumption of the screening properties in the charge chan- 
nel of the particles participating in the plasma analogy 
based on ^'m state. As the distance is increased there 
are less of them and the breakdown of the description in 
terms of dipoles of neutral fermions at smaller distances 
becomes a possibility. We expect that due to impuri- 
ties there will be patches (islands) of dissociated neutral 
fermions 1241. 
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IV. QUANTUM FLUCTUATIONS AND 
QUANTUM DISORDERING 

A. Introduction 

The two paradigms - models of superfluid disordering 
as applied to our 2+1 dimensional system mean that the 
time evolution is such that (1) meron -antimeron pairs 
are locked on impurities or (2) created and annihilated 
at some later time and therefore making a loop in time. 
The loops in time signify the presence of quantum dis- 
ordering [TJ. We will discuss and detect the presence 
of quantum disordering in the WFs of class (c) in Fig. 
1 by examining how they relate to and incorporate or- 
dinary (not quantum disordering that involves merons 
- vortices) quantum fluctuation phonon contribution in 
this case [9, !§] . We will find that the WFs of class (a) in 
Fig. 1 do not have this contribution. 



B. Quantum fluctuations due to phonons and 
quantum disordering 

The usual [l^ Chern-Simons (CS) field theory ap- 
proach [sl in the random phase approximation (RPA) to 
the bilayer problem at v ~ 1 (which in the neutral chan- 
nel reduces just to the problem of ordinary superfluid 
with only phonon description and contribution) finds the 
following correction to the ^'m state: 



PH 



where p^. 



exp 



pI 



{-I 



E 

k 



PE 



-Pk P-k 



(17) 



pI V^k) 



yn(fc)-VTi(fc) V - 22L 
2 ' "^TT — ' 



V^l = ^exp{— fed} i.e. V-[k) is the interaction in the 
neutral channel, Pe = where m is the electron mass 
and p is the uniform total density. In the small d limit 
V- (fc) = ird and we can expand the expression ^' ph as 



PH 



*11 



(E 

k 



(18) 



where c is a positive constant. The terms after the first 
one represent corrections, in the order of importance, to 
the ^'iii ansatz as d increases. 

On the other hand the WFs of Fig. 1(c) are more 
general as they suggest the form of the correction terms of 
wider class than the one used in the expansion (Eq. fTS)) ') 
with only exception that the class demands equal number 
of pI's and p^'s because in writing down the classes of 
Fig. 1 we explicitly distinguished t's from |'s and fixed 
the number of fs and j's. 

We can start comparing and relating the first phonon 
correction i.e. 



E- 



pUpI 



(19) 



to a wave function of two neutral fermions (t from |) i.e. 
density operators as in Eq. ([H]) but with a pairing between 
them as in class Fig. 1(c). 

Without pairing we would have 



1 



(771T - V21) 



pHm^)pHmi) (20) 



which is identical to zero (no correction) as can be found 
out in Appendix A. This is an important result and shows 
that there are no quantum fluctuations in the first family 
of WFs discussed in Section III. Besides this analytical 
proof, our statement is further corroborated by the fact 
that the computer-generated two neutral fermion state 
also does not exist - see Ref. [l6j . 
Therefore we continue by considering 



1 



|?7iT - 



:pHvi])pKv2i) 



(21) 



where a = 1 if we take g(z) — for the pairing func- 
tion or a — 2 ii g{z) — -p- For a — 1 the expression in 
Eq. ([20|) reduces to the form of the first phonon contri- 
bution in the long-distance limit with the ^ singularity 
(see Appendix A) and for a = 2 this singularity soft- 
ens to ~ — Inj/cZs} where Ib is the magnetic length (see 
Appendix A). We will consider only these most weakly 
pairing cases; the case 17(2;) = - does not produce correc- 
tion as can be seen in Appendix A. 

Next we consider more than two density operator con- 
structions i.e. more than two neutral fermions construc- 
tions as in Eq. ([5]) but instead of the two decoupled Fermi 
seas we have a pairing between neutral fermions: 



,2 Ilk<lim - VldIlp<qiVpl - Vql) 

Det{g {iTf -TJi)} x 
pT (?7i|) •• •p-^(ry„ J*iii(z|, 
1 



d'^Vnl 



Det{- 



Det{ 



-} 



I VI -VI 



} 



pT(r,iT)---pH??«i)*iii(^T>^iX22) 

where in the second expression we used the Cauchy de- 
terminant identity i.e. 

Uk<liVk^ - VI]) ]\p<q{Vpi -Vql) _ ^^^^ 1 



U^.AV^] -Vjl) 



' V] - VI 



-} 

(23) 

and substituted the pairing function that has lead us 
to the first phonon correction for two paired neutral 
fermions. Immediately we can see that the diagonal 
terms in which pairs of the two determinants are the same 
would make further phonon-like corrections i.e. their su- 
perposition with appropriate coefficients would lead to 



k 



(24) 
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The other non-diagonal terms would lead to more com- 
plicated constructions of four and more neutral fermions 
that should participate in the description of quantum 
disordering i.e. describe the physics beyond phonon con- 
tribution (|24p . Although in some sense we arc talking 
just about a class (a pool) of wave functions that should 
describe quantum disordering we can fix general form, at 
least for small d, of the superposition that should com- 
pletely model the ground state at fixed d 



\l/r 



n=0,2,... 



(25) 



In the long distance limit ((25)) should tend to ((24l) . In 
other words non-diagonal terms in Eq. (|22p should be sub- 
leading to the leading behavior in ([24]) . That this is 
true from the physical point of view we expect that it 
is enough to prove the subleading behavior in the case 
of four neutral fermions (5'|) and that can be found 
in Appendix B. The proof is based on the smallness of 
higher order terms that may appear inside the brackets 
in Ea. p4)) . This is assumed in the RPA approach and 
expected in the small d limit. 

Therefore the quantum Hall physics besides g{z) — 
y/j^ pairing possibility brings or allows the possibility 
oi g{z) = 4r pairing that introduces non-trivial quantum 
corrections i.e. brings another kind of quantum disor- 
dering. The g{z) = ^f-^ accommodates the usual (on 
the level of RPA) superfluid description in which we may 
expect that the disordered phase will break translation 
symmetry. Indeed, the bosonic CS field theories that 
are not based on quantum Hall WFs give this scenario 
of the disordered phase as a charge density wave 'S]. It 
seems, therefore, there are two possible scenarios for su- 
perfluid disordering not in the BKT class for the bilayer 
in the translation symmetry invariant case (without im- 
purities). In the following we will discuss the second pos- 
sibility with g(z) = 4r kind of pairing. 



C. Weak pairing g{z) ~ case and conformal field 
tlieory considerations 

We expect, if the translational symmetry of the ground 
state remains unbroken, that also in the case of pair- 
ing g{z) — the translatory invariant system smoothly 
evolves with the increase of d into the class of wave func- 
tions in Fig.l (d). We would like to know more about 
this class - whether it represents a distinct phase. If we 
take the choice g{z) = -pr and examine the final form of 
the state of Fig. 1(d) when there are no CBs, we are lead 
to its following forms. 



*2 = 1,. }n(^»T-^jT)'n(^fci 



^Jl t<j 



k<l 



= Vet{^^—r}Vet{ iM-ni, (26) 



where to get the last line we used the Cauchy determi- 
nant identity. The neutral part of 4'2 (not carrying a 
net flux through the system as ^'m does) that consists 
of the two determinants can be viewed as a correlator of 
vertex operators of a single nonchiral bosonic field. Ac- 
cording to |25| conformal field theory (CFT) correlators 
not only describe quantum Hall system WFs but also can 
be used to find out about excitation spectrum and con- 
nect to its edge and bulk theories. In this way motivated 
neutral excitations are vertex operators that correspond 
to single- valued WF expressions that multiply ^"2: 



exp{i5i(/)(u', w*)} 



\2Si 



|2<Si 



(27) 



*\S2 



exp{i(52^(u', w*)} 



n.( 



,*\S2 



■a 



Zk] - zii - 



(28) 

where (f)(w,w*) = 4'{w) + (f){w*), 9{w,w*) — 4i{w) — 
(f>{'w*), and (j){w) and (t>{w*) are holomorphic and anti- 
holomorfic parts of the bosonic field respectively. 62 must 
be i because of the requirement of single valuedness. For 
detailed explanations of the bosonic CFT analogies see 
Appendix C. 

If the low-lying spectrum were consisting only of Si = 
^ and (^2 = 5 quasiparticle excitations our system would 
be described by so-called BF Chern-Simons theory or the 
theory of 2D superconductor [2^. The mutual statis- 
tics of quasiparticles - quasiparticles and vortices, in this 
theory is semionic (due to the fact that vortices carry 
half-flux ((2^) quantum) and that this is also the case 
with our excitations can be easily checked via CFT cor- 
relators - see Appendix C. Combining the analysis with 
the charge part (^'m) in which only charge 1 excitations 
are allowed (half-flux quantum excitations are strongly 
confined [l^l) we may come to the conclusion that the 
degeneracy of the system GSs on the torus must be 4 
|26. . ^8] But the expression for the first kind of excita- 
tions (Eq. (l77)) ') allows a real continuum for the value of 
Si exponent including 61 = and therefore we expect 
a compressible (gapless) behavior of the system despite 
the incompressibility of the charge channel and seemingly 
topological phase behavior in the neutral sector. Never- 
theless we expect that in our case BF CS theory is a part 
of the description of the pairing phase in a Lagrangian 
in which there is a quadratic non-derivative term in one 
of the two gauge fields; this allows a branch of gapless 
excitations - see Appendix C for details. 

The question may come why we did not do an analysis 
with the projection to the LLL. Certainly the analysis 
is more involved where "reversed flux part" i.e. complex 
conjugated determinant becomes an operator that acts 
on the rest of wave function. Nevertheless, an analysis of 
the edge excitation spectrum [29[ suggests that it can not 
conform to any description of simple, free CFT theories 
i.e. can not belong to a totally incompressible class, and 
it is very likely that the system is, as it follows from our 
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unprojected analysis, compressible in the neutral chan- 
nel. Therefore it is very hard to distinguish the physics 
of the states in Fig. 1(d) and Fig. 1(c) in the trans- 
latory invariant system that involve pairing of the type 

While we were finishing the writing a numerical study 
(of homogenous WFs in the translatory invariant case) 
appeared 30] that agrees with and complements our con- 
clusions on the nature of pairing. 



V. CONCLUSIONS 

In conclusion, we presented two families of wave func- 
tions that describe two possible ways of homogeneous 
disordering of the quantum Hall superfluid with their de- 
tailed description on the basis of the dual (quasiparticle) 
picture of the quantum Hall effect. We also presented 
detailed analysis of the disordering in the translation in- 
variant system on the basis of insights into the pairing 
function of quasiparticles - neutral fermions. A class of 
candidate wave functions was clearly connected with the 
formalism that we find in other (Chern-Simons) theories, 
and the pairing function g{z) ~ 4r was extracted as a 
clear choice that incorporates quantum disordering and 
that will describe the system if it does not transform into 
a CDW (charge density wave) inhomogeneous solution. 
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APPENDIX A 



yields 



We want to prove 



d mi d r]2i 



(»7iT - V2l) 



pHmi)pHv2i) ^ (Ai) 



After switching to the Fourier space, p'^i'i]) 
explifcrf}, the l.h.s. becomes 

27r^y" d^ri^exp{ikrf}plp^_^ 



(A2) 



The angle part of the integration with the help of the 
table integral [3l|: 

/ exp{j/3cosa;} cos{na;} = i"7rJ„(/3) (A3) 



d rj — eKp{ikff} — / dr\i'Kj\(kr) — mJ\{—kr)\ 
V Jo 

I /"°° J rdJo{kr) dJo{~kr)- 



2n 



(A4) 



where we used notation |fc| = fc and in the last line the 
identity for the Bessel functions: Jq — —Ji- On the other 
hand 

2 1 7^ „ 27r 

d rj— exp{—ikri} — —i — , (A5) 
rj k 



and therefore Ea. (|A2p can be written as 



^Y^l d'^V-(^^P{ikr]}plp^k+ I d^V-exp{-ikf]}pl^pi 



^'^^'T.wMP-k - P-kPi) = QED. 
k I'^l 

(A6) 



Next we want to evaluate 



d'-ni] / d'mi-. T^p^{m\)p^{m) (A7) 

Again this reduces in the Fourier space to 

2^Y.j ^'^^ eM^rj}py-u (A8) 

In the case of a = 1 as usual for the real Coulomb inter- 
action in 2D the integral is 



1 



27r 



d'^rj— exp{ikr\ = 2-k / drJ^ikr) = — (A9) 
m Jo k 



In the case of a = 2 we have 



1 



d rj- — exp{ifcr} — 2t: I dr 



Mkr) 



(AlO) 



The integral needs a cut-off at small distances (otherwise 
diverges) which should be included in our effective de- 
scription and as usual can be taken to be ?b (magnetic 
length distance). Therefore, instead of Eq. (|A10p we have 

2n r d/-^^^ = 27rKo{lBk) (All) 
Jo + Ib 

In the small momentum limit we can approximate: 

/^o(^)«-ln{^} + o(z) (A12) 

and therefore our first phonon-like correction in this case 
of pairing is 



Y,HHkiB}py 



(A13) 
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For the case of pairing g{z) = i we have 



(?7iT - miY 



P^MpKmi) (A14) 



which reduces to the solving of the following Fourier 
transform 



1 



77- 



(A15) 



With the help of Eq. (|A3p we have for the value of the 
integral: 



1 



dr—{—TTj2{kr) — 7rJ2(— /cr)} 



We may use the table integral [3l| 



(A16) 



for — 1 < Req < Rev — 5 to find out that the value of the 
integral does not depend on k i.e. 

/ d^T]-^ cxpUkf} = -vr (A18) 
J V 

Therefore the phonon-like correction in this case is pro- 
portional to 

k 

and in the real (coordinate) space this becomes: 
d^rip^T])p^{T])^ III = 

E / d^rjS^rt- z,^)S\r^~ Zj^i)^!/ Ill ^ 

Y^S^{z,^~z, ^1)^111^0 (A20) 
hi 

i.e. no correction at all. 



APPENDIX B 

We consider non-diagonal (non-phonon-like) correc- 
tions that come from the description of quantum disor- 
dering by the class of WFs in Fig. 1(c) when the pairing 
is fixed to be g{z) = -y/^, i-C. non-diagonal terms of 
Eg. ip^ with n = 4. We want to prove the subleading 
behavior with respect to the diagonal terms as the one 
with 7711 = ?7i,773| = 773,?72i = and 7741 = m in 



d m d r]3 d 'i]2 d -qi 



1 



1 



\m - ml \v3 - Vil 
pHm)pHm)p^{m)pKm) 



k 

k 



(Bl) 



of the following non-diagonal term: 



(?7i - m) im - Vi) 

/ {m - m) iv3 -m) 
V W - ^4*) (%* - %*) 

pHm)pHm)pKm)p^{m)- 

(B2) 



The non-diagonal terms by their forms should describe 
different processes from the phonon contributions i.e. 
from those as {p\ pi,. ) • • • (pl^ pit„ ) for arbitrary fc's. 

1 ,1 ^ n. ^ n 

In the long-distance approximation we will argue that 
the non-diagonal term fEa. (|B2p ) carry less importance 
that the phonon contribution with the same number of 
density operators. 

Introducing 77 = 771 - 7/4, 77 773 - 772, 77- ^Vi-Vi^ and 
77-|_ = 771 + 773 we can rewrite Eq. (jB2l) as 



J2 Jd^rjJ d^v J d'v^ J d' 



{11- +fj) [ri-ri-) 



77 77 



— exp{i?7lfci} exp{i77^fc3} 

X exp{i(i7/+ - i7f_ - fi)k} 
X exp{i(i7/+ + ^77^ - 7/)fc} 



T T i i 
^^fci Pk, PI Pk 



(B3) 



The 7/+ integration brings the constraint k + k + ki + k^ — 
0. Then the remaining 77- integration gives the following 
contribution, 



1 1 r7_ - - ~ - 

7 —37 7 7 exp{i— (fci - fc3 - fc + fc)} 

(77_ -^77) (77-7;_) 2 



27r 1 ^ ^ :^ ^ ;^ 

—i— — r[exp{77)(fc3 -I- fc)} — exp{— 777(fc3 + k)}] 

\k + kj,\ri + V 



(B4) 



where we used the constraint. Therefore the contribution 
is proportional to 



1 



77 rj 



y I d^v I d^v, . 

Ik + kslJ J iv + v) V 

k3,k,k ' 

[exp{i77(fc3 + k)} — exp{— 77/(^3 + k)}] 
X exp{— 777fc} exp{— 777fc} 

"^Plk-k-k. pI pi pi (B5) 

In the long-distance limit | fc + fc3 1 but that does not 
cancel the part of the 2D volume in the integration mea- 
sure like in the phonon contribution (that would damp 
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the contribution) but is canceled by the difference of the 
exponentials in the same limit in Eq. (|B5p . There is only 
one more factor i.e. -^r^ that can bring: the momen- 

r]+ri o 

turn inverse contribution but this only enforces k k 
i.e. {plp^f.)"^ without a significant coefficient. This will 
only give the next order contribution inside the brackets 
in Eg. ip^ which for small d, and as usual in the RPA 
approach, we can neglect. 



APPENDIX C 

We will give a more general view of the CFT analo- 
gies of so-called [32] doubled CS theories to which BF 
CS theory belongs. In the work of Freedman et al. 
f[3^) BF CS theory was classified as the low-energy the- 
ory of the deconfined phase of Z2 gauge theory. There 
also SU{2)i X SU{2)-^ doubled CS theory was consid- 
ered. For the detailed description of these theories the 
reader should consult Rcfs. |32 andU^. Here we will, by 
writing down relevant CFT correlators, demonstrate the 
analogies between non-chiral - complete CFTs and these 
doubled CS theories. 

First we will consider SU {2)i x SU{2)^ case. The 
possible wave function with coordinates of two species 
. . . , ZATj, for which there are equal number of | 's and 



i's: iV| 



Ni and Nt^ + = N, is 



(CI) 



We use the following correlator of vertex operators of a 
bosonic field S: 



< exp{i(3(j){zi, zl)} exp{-il3(j){z2, Z2)} >= 



1 



\Zl - Z2 



12/32 



(C2) 



If a = ^ we can rewrite our wave function as 
V 2 

^' =< exp{ia(j>{zi, z^)} exjp{ia(l){z2, ^2)} • • • 

exp{-iQ;<^(zjv, ^at)} > (C3) 

and define: 

0(z,z*) =0(z)+. 



e{z,z*)^4>{z) 



(C4) 
(C5) 



Inserting a neutral pair {wi and ^2) of exp{z(5i(/)(u;, w*)} 
vertex operators or ex.-p{i52<t>{w , vertex operators we 
can conclude that these insertions correspond to multi- 
plying the wave function (Eq. (ICip ) by: 



exp{i(5i(/)(z/;, w*)} 



\2s1a 



(C6) 



exp{i(52^?(w, w*)} 



w 



*^(52 



" Y{^{z^l-W) 



(C7) 

(The general formula for the many vertex correlator can 
be found, for example, in Ref. [H.) The single- valuedness 
of the WFs demands ^2 = If we take also = = 
52 — 5 then 

< exp{i(50(wi, wJl)} exp{— i(5(/)(u'2, W2)} ^ 
eicp{i56(wz, w'^)} exp{— i(56'(w4, w\)} > = 
1 1 



\wi - 


W2\''' 


\W3 - 




{wi - 




(W2 - 




K- 




K- 


- wiY' 




-wlf 






{wi 


- WiY' 


{W2 


- 



(C8) 



and the mutual statistics between any of two particles of 
different kinds ((13), (14), (23), or (24)) is fermionic. To 
see that, for example, for the (13) pair we send 2 towards 
4 and switch wi and coordinates. 
In our case of the quantum Hall bilayer. 



nfc<i l^feT - ^iT P Up<q \Zpl - Zqi I 



Det{- 



-}Det{- 



1 



' Zk] - zii - 



-}, 



(C9) 



The same analysis as above will fix a = 1 and (5 = ^ so 
that in this case the mutual statistics is semionic just as 
it should be in the BF CS field theory. 

The BF CS theory of a 2D superconductor is [2^ 



(CIO) 



where and = 0,1,2 are gauge fields; the first 

term is the CS term and j'^ and j^^^J. — 0,1,2 repre- 
sent quasiparticle and vortex density-currents. The La- 
grangian encodes in the i coefficient mutual semionic 
statistics between the two excitations in 2D superconduc- 
tor - any time quasiparticle encircles vortex it gets the 
Bohm - Aharonov phase tt, because vortex corresponds 
to the half-flux quantum excitation in the paired system. 
Higher order in derivatives i.e. Maxwell terms ~ (9 x a)^ 
and ~ {d X b)^ are present in the description of the or- 
dinary (s - wave) - gapped 2D superconductor and can 
describe the plasmon modes that are gapped - see Ref. 
26. In our case, because from CFT analogies fEq. (jC6p 
and Eq. (jC7P ) we find that 5i can be continuous and cor- 
respond to a branch of gapless excitations, we expect a 
quadratic in one of the gauge fields, without derivatives, 
term to describe such a behavior. For example, if we add 
a term quadratic in b {b^b^ with the {d x a)^ Maxwell 
term present) our classical equations of motion will be : 
Aa — and d x b — 0. They describe gapless behav- 
ior (Goldstone mode) in one gauge field and associated 
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quasiparticle description, and incompressible behavior in and we see explicitly mutual fermionic statistics.) 
the other. 

(The SU{2)i X SU{2)^ theory can be described by the 
following Lagrangian 

^ef'^Xd^ax ~ a^r - (Cll) 
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